Resonant two-photon single ionization of two atoms 
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Resonant two-photon ionization in a system consisting of two spatially well-separated atoms is 
studied. Due to two-center electron-electron correlations, the ionization may also proceed through 
photo-excitation of both atoms with subsequent interatomic Coulombic decay. We show that this 
channel may dominate the photoionization process and qualitatively change its dependence on the 
field intensity and the spectra of emitted electrons. 
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The development of tunable laser sources in the 1970s 
stimulated in-depth studies of the behavior of atoms 
and molecules exposed to intense resonant electromag- 
netic fields P, 0]. While the emission of fluorescence 
light can provide detailed atomic structure information, 
resonance-enhanced multiphoton ionization has become 
a valuable tool for probing dynamical aspects of atoms, 
clusters, biomolecules and chemical reactions 

Resonant two-photon ionization of atomic hydrogen 
has been examined thoroughly by theoreticians as a 
prime example which allows a detailed understanding of 
the underlying field-induced dynamics. Due to the reso- 
nantly driven bound-bound transition, the photoelectron 
spectra display Autler-Townes doublets and the total ion- 
ization probability exhibits a step-wise temporal evolu- 
tion Other nonperturbative effects in strong laser 
fields were studied as well (see e.g. [3, S] and references 
therein). Corresponding experimental studies are becom- 
ing feasible nowadays due to the advent of brilliant high- 
frequency photon sources such as advanced synchrotron 
beam lines and free-electron lasers. 

Ionization may also occur through resonant photoex- 
citation of an autoionizing atomic level, which subse- 
quently decays via Auger emission. In recent years this 
type of ionization mechanism, which relies on electron- 
electron correlations, has been investigated intensively, 
both theoretically [6|-l8[ and experimentally, in systems 
consisting of two (or more) atoms, such as noble ga s 
dimers [§, , metal oxides [ll| or water molecules . 
Here the radiationless deexcitation of one atom is ac- 
companied with the ionization of another neighboring 
atom. This process is commonly referred to as inter- 
atomic Coulombic decay, ICD 

In this Letter, we study resonant two-photon ioniza- 
tion of a two-atom system which is exposed to an exter- 
nal electromagnetic field. In this situation, each atom 
is subject not only to the influence of the field but also 
to the interaction with the neighbor atom. Hence there 
are two pathways for ionization: either directly via reso- 
nant two-photon absorption at a single center, or medi- 
ated by ICD when each atom has been excited by single- 
photon absorption (see figure [1]). We shall show that the 
second mechanism can be remarkably efficient and may 



even dominate the direct channel by orders of magnitude. 
Moreover, the two-center channel exhibits very peculiar 
features with respect to its field dependence, temporal 
development and photoelectron spectra. In contrast to 
the two-center resonant photoionization considered in , 
the present process occurs for two identical atoms. It can, 
in principle, also exist in the case of two different atomic 
species provided that the latter ones possess a common 
dipole-allowed transition frequency. 

Let us consider a system consisting of two atoms (A 
and B) separated by a sufficiently large distance R such 
that their individuality is basically preserved. Let each 
atom have an excited state reachable from its ground 
state by a dipole-allowed transition and the energy be- 
tween these states be the same for both atoms. The 
atoms, which are initially in their ground states, are em- 
bedded in a resonant electromagnetic field. As seen from 
figure [U already in the absence of the neighbor, each 
of the atoms can be ionized by absorbing two photons. 
However, due to ICD, photo ionization in the system of 
two atoms acquires interesting and qualitatively new fea- 
tures. 

Assuming the atoms to be at rest, we take the position 
of the nucleus of atom A as the origin and denote the 
coordinates of the nucleus of atom B, the electron of 
atom A and that of atom B by R, ri and r2 = R + £ , 
respectively, where ^ is the position of the electron of 
atom B with respect to its nucleus. 

The total Hamiltonian describing two atoms in the ex- 
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FIG. 1: Scheme of resonant two-photon single ionization of two 
atoms. In order not to overload the picture, ionization channels are 
shown for atom A only. Ionization of atom B proceeds analogously. 
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ternal electromagnetic field reads 

H = Ho + Vab+ W, (1) 

where Hq is the sum of the Hamiltonians for the nonin- 
teracting atoms A and B, Vab the interaction between 
the atoms and W = Wa + Wb the interaction of the 
atoms with the electromagnetic field. 

For electrons undergoing electric dipole transitions the 
interatomic interaction Vab reads 

i?3 ^ ' R 

Here, a = J'li^jl^ - 'iRiRj/R^), (3 = rii£.j{Sij - 
RiRjl R^), rii and S^j G {x, y, z}) are the coordinates 
of the electrons and a summation over the repeated in- 
dices is implied. Further, k = ujfi/c where ujfi is the 
transition frequency and c the speed of light. Note that 
atomic units (a.u.) are used throughout unless otherwise 
stated. 

The electromagnetic field will be treated as a classical, 
linearly polarized field, described by the vector potential 
A(r, t) — cFq/loq cos (coot — ko • r), where luq — cko and 
ko are the angular frequency and wave vector, and Fq is 
the field strength. The interaction W then reads 



(3) 



where Pj is the momentum operator for the j-th electron. 

In what follows we shall assume that the electromag- 
netic field is not too strong such that field-induced tran- 
sitions from bound states to the continuum are weak and, 
besides, one can neglect the infiuence of the field on the 
continuum. This also implies that only one of the cen- 
ters, either A or B, will be ionized with a non-negligible 
probability. 

However, even a relatively weak, but resonant, electro- 
magnetic field can very effectively couple bound states 
belonging to the same center. Therefore, for each cen- 
ter we first construct field-dressed bound states. Taking, 
as an example, the center A, assuming that the field is 
switched on suddenly at U — 0, and using the rotating- 
wave approximation we obtain (see also P, 

ip+{t) = --^^-^—j{{z^+ujo-ei)exp{-izft) 
- (z^ + Wo - £i) exp(-iz^t)) ^0 
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z^ — z: 
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— (exp(-zz^t) - exp(-iz^t)) <^o 



((z^ - £o) c-xpi-iz'^t) 



(z_ — Eq) exp(— e'xjp{—iujat)ipi. (5) 



In these expressions (/Sq with an energy Eq and ipi with an 
energy ei are the ground and excited states, respectively, 
of center A. Further, 



ei - Wo ± riflO : 



(6) 



where = y (ei — eo — wo)^ + 4 | Woi | is the Rabi fre- 
quency and Wij = (ifi |Fo ■ pi/(2a;o)| ifj) G {0, 1}). 
The states @ and JS]) are orthogonal to each other and to 
the continuum states of center A, and are normal- 

ized to unity. At t = these field-dressed states reduce 
to (^0 and ipi, respectively. In the above description we 
have neglected the spontaneous radiative decay of the 
excited state i^i which in our case is justified as long as 
|W^oi| ^ Tr, where F,. is the radiative width of 

Expressions for the field-dressed bound states x±(t) 
on the center B with the corresponding quasi-energies 
can be obtained from (|4|)-([6]) by straightforward re- 
placements. 

Our consideration of the photoionization process will 
be based on the S'-matrix formalism. Concentrating for 
the moment on the description of ionization of atom A 
we write down the transition matrix element 



dt 



Vab + Wa 



(7) 



Here 5" is a solution of the Schrodinger equation for the 
total Hamiltonian H and is a final channel corre- 
sponding to ionization of A. Note that neither Wa nor 
Vab are included in the Schrodinger equation for tp^. 

Taking into account that (i) the electromagnetic field- 
induced transitions to the continuum are weak, (ii) the 
field does not distort continuum states, (iii) the inter- 
atomic interaction at large distances is weak, and (iv) 
initially (at U — 0) both atoms are in the ground states, 
one can approximate '^{t) by ip+{t)x+{t)- Besides, the 
final state tpf {t) can be represented by either one of the 
states LppX+{t) and ippX-{t). Here ipp is the continuum 
state of center A with an asymptotic momentum p and 
energy 

As a result, ionization of center A is described by the 
following transition amplitudes 



Sp,+ {t) = -'^ J dt\ippx+ 



Vab + Wa 



Sp.-{t) = -i / dt' 



Vab 



(8) 



Since the final states ippX+{t) and ^pX-{i) are orthogo- 
nal to each other, these amplitudes add up incoherently 
and for the probability of ionization of center A we obtain 



PA{t) = j d^p{\Sp,+{t) \^ + \Sp 



-it) 



(9) 
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Similar expressions hold for ionization of center B. Note 
that the time integrals in ^ are easily taken analytically. 
However, the resulting expressions are somewhat lengthy 
and will be given elsewhere. 

Let us now turn to the discussion of some results fol- 
lowing from the above expressions. Perhaps one of the 
most appropriate objects, where the process under con- 
sideration may occur, is a helium dimer exposed to a 
resonant field. However, for simplicity we restrict our at- 
tention here to the most fundamental two-atomic system 
consisting of two hydrogen atoms, which are initially in 
the ground state and are irradiated by an electromag- 
netic field with frequency resonant to the atomic ls-2p 
transition. The internuclear distance R is assumed to be 
such that the changes in the atomic levels, caused by the 
(full) interatomic interaction, remain smaller than even 
the natural width Tr of the 2p level in atomic hydrogen, 
which is fulfilled at i? - 20 a.u.. Since we also suppose 
that |Woi| ^ (which in case of hydrogen holds for 
Fq ~ 10~^ a.u.) these changes are negligible compared 
to those induced by the electromagnetic field. Therefore, 
one can indeed treat the ionization process by regard- 
ing the two-atomic system as consisting of two individ- 
ual hydrogens whose atomic properties remain basically 
unchanged. 

For not too strong electromagnetic fields and not too 
large separations between the atoms the two-center chan- 
nel can strongly outperform the direct channel in photo 
ionization. Indeed, neglecting the retardation effects, 
one can show that for ionization of center A the rela- 
tive strength of the former with respect to the latter is 
determined by the ratio (^dB/{R^Fo)^ , where ds is the 
dipole moment of the bound-bound transition in center 
B. The case, when this ratio substantially exceeds unity, 
is illustrated in figure [2] where we present the probabil- 
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FIG. 2: Probability for single ionization of two hydrogen atoms 
separated by a distance R = 25 a.u. given as a function of the 
electromagnetic pulse duration for zero detuning. The interatomic 
vector R is directed along the field polarization Fq. Dot, dash and 
solid curves show the results for Fq = 2 X 10~^, 5 X 10~^ and 10~^ 
a.u., respectively, a) the contribution of the direct channel only, b) 
the contribution of the two-center channel. 



ity PA{t) + Psit) for single ionization of two centers as 
a function of time. Since electron transitions into the 
continuum states occur only from the excited states the 
probability shows a non-monotonous behavior in which 
time intervals, when the ionization probability rapidly in- 
creases, are separated by intervals, when the probability 
remains practically constant, refiecting oscillations with 
the Rabi frequency of the electron populations between 
the ground and excited states in a resonant electromag- 
netic field. 

Such a staircase behavior, related to the Rabi oscilla- 
tions, is inherent for both the direct [J and two-center 
channels of ionization (see figure[2^ and[2j), respectively). 
However, since in the latter case for ionization to occur 
both centers have to be in the excited states, the "stairs" 
in the time development of the probability for ionization 
via ICD are more pronounced. 

Note also that the probability for ionization via ICD 
may demonstrate a behavior which appears counterin- 
tuitive at first glance: the weakest field can lead to the 
highest instantaneous value of this probability (see fig- 
ure Wp, pulse durations in the intervals around 75 and 
175 ps). However, a simple analysis shows that such a 
behavior is in fact the consequence of (i) the dependence 
of the Rabi frequency on the field intensity and (ii) the 
field-independence of the probability for ionization via 
the two-center channel, temporally averaged over the in- 
verse of the Rabi frequency. 

The direct and two-center channels of photo ioniza- 
tion are in general characterized by different dependences 
on the field intensity. For instance, in the case of ion- 
ization by a monochromatic field, which is exactly res- 
onant to the bound-bound transitions, the population 
probabilities for the excited bound states (averaged over 
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FIG. 3: Probability for single ionization of two hydrogen atoms 
separated by a distance R = 25 a.u. given as a function of the 
electromagnetic field strength for a pulse duration of 100 ps. R 
is directed along Fq. Solid curve: ionization by a monochro- 
matic field with zero detuning. Dash curve: ionization by a non- 
monochromatic field with bandwidth of 20 meV and central fre- 
quency equal to the ls-2p transition frequency in hydrogen. 
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the inverse of the Rabi frequency) do not depend on 
the field intensity. As a result, the averaged probabil- 
ity for photo ionization via the direct channel, which in- 
volves absorption of an additional photon, is proportional 
to the field intensity. In contrast, the averaged prob- 
ability for ionization due to the two-center mechanism 
is intensity-independent because it relies only on reso- 
nant bound-bound transitions and two-center electron- 
electron correlations. These two limiting dependences 
and a smooth transition between them are illustrated by 
the solid curve in figure El Similar results hold also for 
non- monochromatic fields (see dash curve in figure [3|) . 
At higher intensities, where the direct channel is more 
efficient, the ionization probability demonstrates a lin- 
ear growth with intensity whereas at smaller intensities, 
where the two-center mechanism dominates, the ioniza- 
tion probability becomes almost a constant. The wiggles 
seen in figure [3] result from the Rabi- flopping dynamics. 

In a resonant electromagnetic field each of the ground 
and excited levels of the centers A and B split into 
two sub-levels. Therefore, the direct ionization channel 
leads in general to two Autler-Townes lines (per atomic 
species) in the energy spectrum of the emitted electrons 
located at Sp ^ z±. + 2ujo (see figure H^) . Because of 
the splitting electron transitions between dressed bound 
states are characterized by three different transition ener- 
gies. When these energies are transferred to the neighbor 
field-dressed center via ICD, one would obtain in gen- 
eral six energy lines (per atomic species) for the emit- 
ted electrons. For the emission from center A they are 
given by Sp — z±. + 2ujo, Sp = z±. + — + 2cjo and 
Ep = z:^ + zf — z-^ + 2ujQ. However, for identical centers 
only four of them have different energies (see figure |3]d). 

Since in the transition amplitude >5'p,+ both direct and 
ICD ionization channels may lead to the same final states, 
they can interfere. As was mentioned, for identical atoms 
the direct and ICD channels result in two and four emis- 
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FIG. 4: Energy spectrum of electrons emitted in the process of 
single ionization of two hydrogen atoms separated by a distance 
R = 25 a.u. at Fo = 10~* a.u., zero detuning, pulse duration of 10 
ps and R directed along Fq . Panels a) and b) show the partial con- 
tributions of the direct and two-center mechanisms, respectively. 
Panel c) displays the total spectrum including the interference of 
the two channels. 



sion lines, respectively, in the energy spectrum of the 
emitted electrons. The former two coincide with the two 
central lines from the latter four. Therefore, for these two 
lines interference may occur. Indeed, comparing figures 
|3^,|Dd andSt we see that the total contribution of the two 
channels is not simply equal to the sum of their partial 
contributions. In particular, destructive and constructive 
interferences of these contributions in the amplitude Sp^+ 
occur, respectively, for the left and right central peaks. 



In conclusion, resonant two-photon single ionization of 
a two-atom system was considered. For simplicity and 
in order to clearly reveal the basic underlying physics re- 
sults for two hydrogen atoms were shown. It was found 
that, in a certain range of interatomic distances and ex- 
ternal field strengths, the ionization is dominated by the 
two-center channel involving ICD. In this case, the mean 
ionization yield - averaged over one Rabi oscillation - 
becomes independent of the applied field intensity. The 
instantaneous ionization probability shows a step- wise in- 
crease with time and a non-monotonous dependence on 
the field strength. The photo-electron spectrum com- 
prises four lines. Its asymmetry reflects the interference 
between the direct and two-center ionization channels. 
An experimental observation of the predicted effects can 
be possible utilizing helium dimers, which were recently 
used successfully for studies of various aspects of ICD 
0. 



[1] P. L. Knight and P. W. Milonni, Phys. Rep. 66, 21 (1980) 
[2] M.V. Fedorov and A.E. Kazakov, Prog. Quant. Electr. 
13, 1 (1989). 

[3] J. C. Miller et al, Phys. Rev. Lett. 45, 114 (1980); G. 

von Helden et al., Phys. Rev. Lett. 79, 5234 (1997); E. 

Nir et al, Chem. Phys. Lett. 355, 59 (2002). 
[4] K. Rzazewski and J. H. Eberly, Phys. Rev. Lett. 47, 408 

(1981); G. S. Agarwal et al, ibid. 48, 1164 (1982). 
[5] M.G. Girju et al, J.Phys. B 40, 4165 (2007) 
[6] L. S. Cederbaum et al., Phys. Rev. Lett. 79, 4778 (1997). 
[7] V. Averbukh, I. B. Miiller, and L. S. Cederbaum, Phys. 

Rev. Lett. 93, 263002 (2004); K. Gokhberg et al, Euro- 

phys. Lett. 72, 228 (2005) 
[8] B. Najjari, A. B. Voitkiv and C. Miiller, Phys. Rev. Lett. 

105, 153002 (2010) 
[9] T. Jahnke et al., Phys. Rev. Lett. 93, 083002 (2004); Y. 

Morishita et al, ibid. 96, 243402 (2006). 
[10] N. Sisourat et al. Nature Phys. 6, 508 (2010); T. Haver- 

meier et al, Phys. Rev. Lett. 104, 133401 (2010). 
[11] A. Kay et al. Science 281, 679 (1998); F.J. Garcia de 

Abajo et al, Phys. Rev. Lett. 82, 4126 (1999) 
[12] T. Jahnke et al, Nature Phys. 6, 139 (2010); M. Mucke 

et al, ibid. 6, 143 (2010). 



